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Abstract: We describe how the coupling of the gluonic Polyakov loop to quarks solves 
different inconsistencies in the standard treatment of chiral quark models at finite temper- 
ature at the one quark loop level. Large gauge invariance is incorporated and an effective 
theory of quarks and Polyakov loops as basic degrees of freedom is generated. From this 
analysis we find a strong suppression of finite temperature effects in hadronic observables 
below the deconfinement phase transition triggered by approximate triality conservation 
in a phase where chiral symmetry is spontaneously broken (Polyakov cooling). We also 
propose a simple phenomenological model to describe the available lattice data for the 
renormalized Polyakov loop in the deconfinement phase. Our analysis shows that non per- 
turbative contributions driven by dimension-2 gluon condensates dominate the behaviour 
of the Polyakov loop in the regime T c < T < 6T C . 
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1. Introduction 

Pure gluodynamics formulated using the imaginary time formalism of finite temperature 
field theory has an extra discrete glogal symmetry Z(iV c ), which is the center of the usual 
gauge group SU(iV c ). A natural order parameter for the transition from the confining phase, 
where Z(iV c ) is preserved, to the deconfining phase, where this symmetry is spontaneously 
broken, is the traced Polyakov loop (for a comprehensive review see e.g. defined by 

L(T) = (tr Q(x)) = (^tr P tf /T dx A (x,x )^ ^ ? ( L1 ) 

where ( ) denotes vacuum expectation value. Aq is the gluon field in the (Euclidean) time 
direction. There have been many efforts in studying the Polyakov loop. A perturbative 
evaluation of the Polyakov loop was carried out long ago ||] at high temperatures. An 
update of this calculations was presented by us in j|] . Different renormalization procedures 
have been proposed on the lattice simulations more recently [Q, ||]. 

In full QCD, i.e. with dynamical fermions, the Polyakov loop appears to be an ap- 
proximate order parameter, as lattice simulations suggest ||. This may look a bit puzzling 
since the center symmetry is largely broken for current quarks. However, as discussed in 
Q in the context of chiral quark models, the relevant scale stemming from the fermion 
determinant is in fact the constituent quark mass, generated by spontaneous chiral sym- 
metry breaking. Thus, one expects large violations of the center symmetry to correlate 
with chiral symmetry restoration. 

We have analyzed the role of large gauge symmetry in a similar framework ||] yielding 
a unique way of coupling the polyakov loop to effective constituent quarks. In Ref. ||, 
we have also proposed a model to describe the available lattice data for the renormalized 



Polyakov loop in terms of the dimension-2 gluon condensate [10, 11]. This model also 
describes consistently the lattice results for the free energy fl! 
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2. Large gauge transformations 



In the Matsubara formalism of quantum field theory at finite temperature the space-time 
manifold becomes a topological cylinder. In principle, only periodic gauge transformacions 
are acceptable since the quark and gluon fields are stable under these transformations: 

g(x,x ) =g(x,x + (2.1) 

where (3 = 1/T. In the Polyakov gauge, where OqAq = 0, Aq is a diagonal and traceless 
N c x N c matrix. We can consider the following periodic gauge transformation 

g (xo) = j 2 ™° A /P , (2.2) 

where A is a color traceless diagonal matrix of integers. Note that it cannot be considered 
to be close to the identity, and in that sense we call it a large gauge transformation. The 
gauge transformation on the Aq component of the gluon field is 

A ^A + jA, (2.3) 

and so gauge invariance manifests as the periodicity of the diagonal amplitudes of Aq of 
period 2ir/ (3. This invariance is manifesly broken in perturbation theory, since a periodic 
function is approximated by a polynomial. Nevertheless, we can implement this large gauge 
symmetry by considering the Polyakov loop or untraced Wilson line as an independent 
degree of freedom, Q(x), which transforms covariantly at x 

n(x)^g- 1 (x)n(x)g(x), (2.4) 

and, in the Polyakov gauge, £l(x) = gtP-AoW f it becomes gauge invariant. 

Fermions break the center symmetry of the gauge group, which is present in all the 
pure gauge theories. That means that we can only consider periodic gauge transformations 
(see Eq. ( |2.1| )). In pure gluodynamics at finite temperature one can smooth this condition, 
and consider aperiodic gauge transformations: 

g(x,x + 0) = zg(x,x ) , z Nc = 1 . (2.5) 

Note that z is not an arbitrary phase but an element of Z(iV c ). An example of such a 
transformation in the Polyakov gauge es given by 

g(x ) = e i2 ™oA/AW3 } (2.6) 

for which z = e l2w ^ Nc . The corresponding gauge transformation of the Aq field and the 
Polyakov loop 0, is 

A ^A + —A, Q^zn (2.7) 

We observe that 0, transforms as the fundamental representation of the Z(A r c ) group. From 
Eq. (|2.7| ) we deduce that (O) = z(Q) and hence (f2) = in the center symmetric or confining 
phase. More generally, in this phase 

{n n ) = for n± mN c , (2.8) 
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with m and arbitrary integer. Obviously, in full QCD the fermion determinant changes 
the selection rule Eq. (J^J^). This violation is large for massless quarks since we expect 
corrections ~ e _M °/ T trf2, and the usefulness of the center symmetry becomes doubtful. 



3. Problems with Chiral Quark Models at finite temperature 

The standard treatment of Chiral quark models at finite temperature presents some in- 
consistencies. To illustrate this point we will use the Nambu-Jona-Lasinio model. In the 
Matsubara formalism we have the standard rule to pass from T = formulas to T / 0, 



/ 



dko 



F(k,k )^iT J2 F$,iuj n ), (3.1) 



where uj n = 2irT(n + 1/2) are the fermionic Matsubara frecuencies. The chiral condensate 
at finite temperature at one loop level with that rule is given by 

/cPk 1 
W ^TFW ' (3 - 2) 

where M is the constituent quark mass. Doing the integral and after a Poisson resumma- 
tion, we have the low temperature behaviour 

7r z ^-^ n 

71=1 

Lo " T fc-fD-ir ^y /2 e -M/T ; (33) 

n=l \ n7r / 

where we have used the asymptotic form of the Bessel function K n (z). This formula can 
be interpreted in terms of the quark propagator in coordinate space 

f d A k e~ ik - x ,._ M 2 K 1 (V-M 2 x 2 ) , n A , 

so that at low temperature we get S{x,i(5) Lo ~ T e~ M ^ T , which represents the exponential 
suppression for a single quark at low T. This means that we can write the quark condensate 
in terms of Boltzmann factors with mass M n = nM. We have the general result that to 
any observables which are color singlets, quark models calculations at finite temperature 
in the one loop aproximation generate all possible quark states, i.e. 

Q T = T=0 + 0ge -M/T + 0gqe ~2M/T + (35) 

Note that while the term O q corresponds to a single quark state, the next term O qq must be 
a qq diquark state, corresponding to a single quark line looping twice around the thermal 
cylinder. It cannot be a qq meson state because at one loop this state comes from the quark 
like going upwards and then downwards in imaginary time, so that the path does not wind 
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around the thermal cylinder and then it is already counted in the zero temperature term 
qT=o From Eq we obtain 

oo 

<3?>t= E (- 1 ) n (9(^o)g(o))U = i ^. (3.6) 

n=— oo 

Under a gauge transformation of the central type we have q(n(3)q(0) — > z~ n q(n(3)q(0). This 
means that Eq. ([O]) is not gauge invariant, and the quark condensate can be decomposed 
as a sum of irreducible representations of a given triality n. 

Another problem comes from comparison with chiral perturbation theory at finite 
temperature. In the chiral limit the leading thermal corrections to the quark condensate 
for Nf = 2 are given by 

384/^ 

which shows that the finite temperature correction is ./V c -suppressed as compared to the 



(gg)lchPT = {qq)r=o ( l - ^ - -^jji + " ' ) > ( 3 - 7 ) 



zero temperature value. This feature contradicts our result of Eq. (3.6) obtained by using 
the standard finite temperature treatment of chiral quark models. The problem is that 
chiral quark model phase transitions located at T c ~ 240MeV are based on such a formula, 
where multiquark states are excited even at very low temperatures. 

4. The Polyakov loop chiral quark model 

We can formally keep track of large gauge invariance by coupling gluons to the model in 
a minimal way. As we said in section ^, a perturbative treatment of the Aq component 
of gluon field manifestly breaks gauge invariance at finite temperature, and we need to 
consider the Polyakov loop as an independent degree of freedom. It appears naturally in 
any finite temperature calculation in the presence of minimally coupled vector fields within 
a derivative or heat kernel expansion j|, 13]. Our approach is similar to that of Q, except 



that there a global Polyakov loop is suggested in analogy with the chemical potential. 
Instead we consider a local Polyakov loop Q(x) coupled to the quarks j|]. From those 



calculations we deduce the rule of Eq. (3.1), but with the modified fermionic Matsubara 
frequencies 

Cj n = 2vrT(n + 1/2 + v) , v = (2vri) _1 log Q , (4.1) 

which are shifted by the logarithm of the Polyakov loop £1 = e t27TU , i.e. v{x) = Aq(x) /(2itT). 
The effect of such a shift over a finite temperature fermionic propagator starting and ending 
at the same point is 

oo 

F(x;x)^ ^2 (-n(x)) n F(x,x + in/3;x,x ) , (4.2) 



instead of the (— l) n factor obtained from the standard rule Eq. (|3.1| ) after using Poisson's 
summation formula and Fourier transformation. 1 To restore gauge invariance we project 
1 This formula can be interpreted saying that in a quark loop at finite temperature, the quarks pick up a 
phase ( — 1) due to Fermi-Dirac statistics, and a non Abelian Aharonov-Bohm factor each time the quarks 
wind once around the compactified thermal cylinder. 
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onto the colour singlet sector by integration of gluons according to the QCD dynamics. 
Effectively this implies an average over the local Polyakov loop with some normalized 
weight cr(ri; x)d£l. Here dQ is the Haar measure of SU(iV c ) and cr(S7; x) the probability 
distribution of Q,{x) in the gauge group. For a general function /(O), meaning a ordinary 
funcion f{z) evaluated at z = O, we have 

K tr - m ) = Ur^ip^'^r^^- <4 - 3) 



where e 1 ^ , j = 1, . . . , N c are the eigenvalues of Q and 

iV c 



•/SU(iV c ) ^ V c r~f 



(4.4) 



By applying the aforementioned rules, the chiral quark model coupled to the Polyakov 
loop corresponds to simply make the replacement 

d -> d - iAo (4.5) 

in the Dirac operator for the quarks, and to consider the partition function 

Z = J DUDQe iTG[n] e ir ^ [U ' n] , (4.6) 

where U is the nonlinearly transforming pion field, DU is the Haar measure of the chiral 
flavour group SU(Nf)n x SU (Nf) l and D£l the Haar measure of the colour group SU (N c ), 
Tq is the effective gluon action whereas Tq stands for the quark effective action. If the 
gluonic measure is left out Ao = and 0, = 1 we recover the original form of the corre- 
sponding chiral quark model, where there exists a one-to-one mapping between the loop 
expansion and the large N c expansion both at zero and finite temperature. Equivalently 
one can make a saddle point approximation and corrections thereof. In the presence of the 
Polyakov loop such a correspondence does not hold, and we will proceed by a quark loop 
expansion, i.e. a saddle point approximation in the bosonic field U, keeping the integration 
on the (constant) Polyakov loop £1. 

By appling these rules to the quark condensate, we deduce in the quenched approxi- 
mation 

oo 1 

(QQ)T= E ^c(-nr)(q(x o )q(0))\ Xo=inP . (4.7) 

n=— oo c 

From Eq. (|2.8| ) we observe that in the confining phase triality is preserved, so that after 



gluon average Eq. (|4.2|) becomes 



oo 

F(x;x)^ Yj ({-^(x)) nNc }F(x,xo + inN c P;x,x ). {4.1 

n= — oo 
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At sufficiently low temperature the distribution of the Polyakov loop becomes just the Haar 
measure, and one can easily deduce the following result 



(tr c (-n)*> suwo 



N c , n = (4.9) 
-1, n = ±N c (4.10) 
0, otherwise (4-11) 



Taking into account this formula in Eq. (4/7), we observe that the inclusion of the Polyakov 
loop not only removes the triality breaking terms, but also the thermal contributions are 
N c suppressed as compared to the zero temperature value, as is expected from ChPT (see 
Eq. ([3,71)), The quark condensate at finite temperature at one loop level and in quenched 
approximation is 

(W)T = (q q )T=o+ 2 -^K 1 (N c M/T) + ... L °^ (gg) T=0+ 4 ^ e "W . (4.12) 

The dots indicate higher gluonic or sea quark effects. Due to the exponential suppression, 
the leading thermal corrections at one quark loop level starts only at temperatures near the 
deconfinement phase transition. We have named this effect Polyakov cooling ||, because it 
is triggered by a group averaging of Polyakov loops. This means that in the quenched ap- 
proximation we do not expect any important finite temperature effect on quark observables 
below the deconfinement transition, and the biggest change should come from pseudoscalar 
loops at low temperatures. This is precisely what one expects from ChPT. 

In order to go beyond the quenched approximation, we will consider the computation 
of the fermion determinant in the presence of a slowly varying Polyakov loop following the 



techniques developped in [13]. Such an approximation makes sense in a confining region 
where there are very strong correlations beween Polyakov loops. The fermion determinant 
can be written as 

Det(i p—M) =e~J rf4 ^' n ) , (4.13) 

where C is the chiral Lagrangian as a function of the Polyakov loop which has been com- 
puted at finite temperature in Ref. |§ in chiral quark models. Using this we can estimate 
the Polyakov loop 2 

J_ (trcg(g) Det(iff-M)) low T 8ttT 2 B _ m/t 

N c (Det(»#-M)) ~ C N^ 6 ' [ ' 

where B is the vacuum energy density, a is the string tension and c is a numerical factor 
which depends on the model. Note that triality is not preserved due to the presence of dy- 
namical quarks, but the relevant scale is the constituent quark mass. So the Polyakov loop 
can be effectively used as an order parameter. In Fig. (|l]) we confront such an exponential 
suppression with unquenched lattice calculations below the phase transition. We observe 
that Eq. ( gig) could be a good approximation below 0.6T C . In any case, lattice data for 



2 The integration can be easily computed by using the formula (tr c f2(:r) tr c f2 1 (y)) = e 



-v\/t 
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Figure 1: Temperature dependence of the renormalized Polyakov loop in units of the critical 
temperature. Lattice data correspond to 2-flavor QCD, and has been taken from p|. The line 
represents our estimation of the Polyakov loop in the low temperature regime, using c = 3 as a 
suitable value for this model-dependent parameter. 



lower temperatures are desirable to do a more precise analysis. For the quark condensate 
we take into account the result of Eq. (4.7), so that 



, (qqDet(i ft-M)) LowT ._ , A ,%ttT 2 B _ 2M/T \ 

{qq)T = (Det(iV-M)) ~ ^ox^l + c'^e-^j. (4.15) 

where, again d depends on the particular model. In other chiral quark models, similar 
results are obtained by replacing 2M — > My (the p meson mass). The Polyakov cooling 
persists although is a bit less effective, and for instance the temperature dependence of the 
low energy constants of the effective chiral Lagrangian becomes Lf — Lf =0 °~ e~ Mv l T . 

Finally, on top of this one must include higher quark loops, or equivalently mesonic 
excitations. They yield exactly the results of ChPT fli~4j| and for massless pions dominate 
at low temperatures. Thus, we see that when suitably coupled to chiral quark models the 
Polyakov loop provides a quite natural explanation of results found long ago on purely 
hadronic grounds. 



5. Polyakov loop above T c 

In this section we focus on the behaviour of the Polyakov loop in the deconfining phase. In 
that phase chiral symmetry is restored and the degrees of freedom are quarks and gluons. 
A perturbative evaluation of the Polyakov loop was carried out in ||] in pure gluodynamics 
to NLO, which corresponds to 0(g 4 ) in the Landau gauge. In the Polyakov gauge Eq. ( |1 . 1| ) 
becomes 

L{T) _ J_ (^(fl/r) _ , _ £±.WA$» + ^-MMM)) +..., (5.1) 
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where we have considered a series expansion in the gluon field. 3 To describe the dynamics 
of the A)( x ) held we use the 3-dimensional reduced effective theory of QCD, obtained from 
the Euclidean QCD action by integrating the non stationary Matsubara gluon modes and 
the quarks || Let Z?oo(k)<5a6 denote the 3-dimensional propagator for the gluon field, 
then 4 

(A 2 ,«> = (iV c 2 - 1)T / — Doo (k). (5.2) 



To lowest order in perturbation theory the propagator becomes DQ (k) = l/(k 2 + m 2 D ), 



where mo is the Debye mass, which to one loop [16 1 writes mo = gT(N c /3 + Nf/6) 



,1/2 



We can compute (Aq) and (Afy by taking derivatives of the vacuum energy density of the 
3-dimensional theory, already computed to four loops in [17]. The contribution from q 2 (Aq) 
starts at 0(g®), while that of g 2 {A 2 ^) starts at 0(g s ). So, the replacement of Eq. (|5.1| ) with 

9 2 {A 2 )- 



L(T) = exp 



AN C T 2 



(5.3) 



becomes correct up to 0(g 5 ). This gaussian ansatz is exact in the large N c limit. We 
obtain 

{ALf = ^ moT _ M£zAfT> ( log - + ») + 0(9 3 } , ( , 4) 



This result can be deduced in two forms: by derivating the vacuum energy density [JTj 
or by identifying Eq. ( |5.3j ) with the perturbative result of Ref. Note that this formula 
holds also in the unquenched theory, since to this order, Nf only appears through the De- 
bye mass. The perturbative contributions to the Polyakov loop at 0(g s ) and C(g 4 ) have 
been displayed in Fig. (||) and compared to the lattice data of the renormalized Polyakov 
loop of Ref. [Q]. It only seems to reproduce these data for the highest temperature value 
6T C . The results of Ref. [0] would provide the 0(g 5 ) and 0{g & ) terms. Unfortunately, 
this perturbative result is obtained in covariant gauges, generating a spurious gauge de- 
pendence beyond 0{g A ). In any case these terms produces a logarithmic dependence with 
temperature so that their contribution is similar to that of the lowest orders. 

It is clear that perturbation theory cannot explain by itself lattice data, and we pro- 
pose to account for non perturbative contributions coming from condensates. We consider 
adding to the propagator new phenomenological pieces driven by positive mass dimension 
parameters: 

D NP(k ) = ^™G . (5.5) 

00 V ' (k 2 +m 2 D ) 2 V 1 

This piece produces a non perturbative contribution to the gluon condensate, namely, 
(A 2 la ) NP = (N 2 — l)Tmg/(87rm£)). If we assume that is temperature independent, 



3 Note that the odd order terms vanish due to the conjugation symmetry of QCD, A^(x) — * —A^(x). 
4 We consider Aq = ^2 a T a Ao, a , with the standard normalization for the Hermitian generators of the 



SXJ(N C ) Lie algebra in the fundamental representation, tr(T a Tb) = <5 a i>/2. 
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Figure 2: Temperature dependence of the renormalized Polyakov loop. Lattice data from 
At the left we plot the prediction of perturbation theory at LO and NLO in pure gluodynamics 
and compare with lattice data from [Q. At the right, the logarithmic dependence of the Polyakov 
loop versus the inverse temperature squared. Fits with a adjustable constant and predicted by NLO 
perturbation theory are displayed. Purely perturbative LO and NLO results for Nj = are shown 
for comparison. 



the condensate will also be temperature independent, modulo radiative corrections. Adding 
the perturbative and non perturbative contributions, we get for the Polyakov loop 

n 2/A2 \P n 2 (A 2 \NP 
- 21 ° gL = ^ATr^ + 2N C T* ' (5 ' 6) 

where (A^ a ) p is given by Eq. ( ^.4[ ). We will rewrite this formula as —2 log L = a+b(T c /T) 2 . 

The lattice data for — 21ogL versus (T c /T) 2 are displayed in Fig. @. The nearly 
straight line behaviour is clear, which means the unequivocal existence of a temperature 
power correction driven by a dimension 2 gluon condensate. If we fit the lattice data by 
using the perturbative value of a up to NLO, i.e. 0(g ), one obtains: 

2.27(6), 2 2 iVP _f (1.00(1) GeV) 2 , N f = 

2.89(8), 9 [ o ' a)T i (0.84(2) GeV) 2 , N f = 2 ' 1 j 

In the second case we use only data above 1.1T C . A fit of the lattice data with a treated 
as a free parameter gives 

-0.24(1), / 1.81(4), 2 2 NP _( (0.89(1) GeV) 2 , N f = 

-0.26(1), 1 2.07(2), 9 { *°> a)T "\ (0.71(2) GeV) 2 , N f = 2 ' 1 j 



Recent analysis of the heavy quark free energy with the model proposed in Eq. (| 
(see [12 1) suggest the possibility that a s at finite temperature has a smoother behaviour 



than the predicted by perturbation theory in the regime T c < T < 6T C . This is in contrast 



with existing analysis 18 1, where authors find a very large value for a s in this regime. 

We can compare our result for the gluon condensate with finite temperature deter- 
minations based on the study of non perturbative contributions to the pressure in pure 



gluodynamics [19|. These results yield for the gluon condensate (0.93 ± 0.07 GeV) in the 



temperature region used in our fit and in Landau gauge, which is in good agreement with 
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Eqs. (5/7) and ( |5.8|) . We also can compare with zero temperature determinations of the 
gluon condensate g 2 (A^ a ) in the Landau gauge and in quenched QCD. From the gluon 
propagator (2.4 ± 0.6 GeV) 2 HjJ and from the quark propagator (2.1 ± 0.1 GeV) 2 fi~l|| . 5 
We observe a remarkable agreement, taking into account that these results refer to different 
temperatures and gauges. 
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